Structural basis for the mechanism of interaction of SARS-CoV-2 B.1.640.2 variant RBD with the host receptors hACE2 and GRP78.

The inflicted chaos instigated by the SARS-CoV-2 (Severe acute respiratory syndrome coronavirus 2) globally continues with the emergence of novel variants. The current global outbreak is aggravated by the manifestation of novel variants, which affect the effectiveness of the vaccine, attachment with hACE2 (human Angiotensin-converting enzyme 2) and immune evasion. Recently, a new variant named University Hospital Institute (IHU) (B.1.640.2) was reported in France in November 2021 and is spreading globally affecting public healthcare. The B.1.640.2 SARS-CoV-2 strain revealed 14 mutations and 9 deletions in spike protein. Thus, it is important to understand how these variations in the spike protein impact the communication with the host. A protein coupling approach along with molecular simulation protocols was used to interpret the variation in the binding of the wild type (WT) and B.1.640.2 variant with hACE2 and Glucose-regulating protein 78 (GRP78) receptors. The initial docking scores revealed a stronger binding of the B.1.640.2-RBD with both the hACE2 and GRP78. To further understand the crucial dynamic changes, we looked at the structural and dynamic characteristics and also explored the variations in the bonding networks between the WT and B.1.640.2-RBD (receptor-binding domain) in association with hACE2 and GRP78, respectively. Our findings revealed that the variant complex demonstrated distinct dynamic properties in contrast to the wild type due to the acquired mutations. Finally, to provide conclusive evidence on the higher binding by the B.1.640.2 variant the TBE was computed for each complex. For the WT with hACE2 the TBE was quantified to be-61.38 ± 0.96 kcal/mol and for B.1.640.2 variant the TBE was estimated to be -70.47 ± 1.00 kcal/mol. For the WT-RBD-GRP78 the TBE -was computed to be 32.32 ± 0.56 kcal/mol and for the B.1.640.2-RBD a TBE of -50.39 ± 0.88 kcal/mol was reported. This show that these mutations are the basis for higher binding and infectivity produced by B.1.640.2 variant and can be targeted for drug designing against it.Communicated by Ramaswamy H. Sarma.

AI-Driven De Novo Design and Molecular Modeling for Discovery of Small-Molecule Compounds as Potential Drug Candidates Targeting SARS-CoV-2 Main Protease.

Over the past three years, significant progress has been made in the development of novel promising drug candidates against COVID-19. However, SARS-CoV-2 mutations resulting in the emergence of new viral strains that can be resistant to the drugs used currently in the clinic necessitate the development of novel potent and broad therapeutic agents targeting different vulnerable spots of the viral proteins. In this study, two deep learning generative models were developed and used in combination with molecular modeling tools for de novo design of small molecule compounds that can inhibit the catalytic activity of SARS-CoV-2 main protease (Mpro), an enzyme critically important for mediating viral replication and transcription. As a result, the seven best scoring compounds that exhibited low values of binding free energy comparable with those calculated for two potent inhibitors of Mpro, via the same computational protocol, were selected as the most probable inhibitors of the enzyme catalytic site. In light of the data obtained, the identified compounds are assumed to present promising scaffolds for the development of new potent and broad-spectrum drugs inhibiting SARS-CoV-2 Mpro, an attractive therapeutic target for anti-COVID-19 agents.

Marine drugs as putative inhibitors against non-structural proteins of SARS-CoV-2: an in silico study.

INTRODUCTION: Coronavirus disease 2019 (COVID-19) is an unprecedented pandemic, threatening human health worldwide. The need to produce novel small-molecule inhibitors against the ongoing pandemic has resulted in the use of drugs such as chloroquine, azithromycin, dexamethasone, favipiravir, ribavirin, remdesivir and azithromycin. Moreover, the reports of the clinical trials of these drugs proved to produce detrimental effects on patients with side effects like nephrotoxicity, retinopathy, cardiotoxicity and cardiomyopathy. Recognizing the need for effective and non-harmful therapeutic candidates to combat COVID-19, we aimed to develop promising drugs against SARS-COV-2. DISCUSSION: In the current investigation, high-throughput virtual screening was performed using the Comprehensive Marine Natural Products Database against five non-structural proteins: Nsp3, Nsp5, Nsp12, Nsp13 and Nsp15. Furthermore, standard precision (SP) docking, extra precision (XP) docking, binding free energy calculation and absorption, distribution, metabolism, excretion and toxicity studies were performed using the SchrÓ§dinger suite. The top-ranked 5 hits obtained by computational studies exhibited to possess a greater binding affinity with the selected non-structural proteins. Amongst the five hits, CMNPD5804, CMNPD20924 and CMNPD1598 hits were utilized to design a novel molecule (D) that has the capability of interacting with all the key residues in the pocket of the selected non-structural proteins. Furthermore, 200 ns of molecular dynamics simulation studies provided insight into the binding modes of D within the catalytic pocket of selected proteins. CONCLUSION: Hence, it is concluded that compound D could be a promising inhibitor against these non-structural proteins. Nevertheless, there is still a need to conduct in vitro and in vivo studies to support our findings.

Molecular insights to the binding interactions of APNS containing HIV-protease inhibitors against SARS-CoV-2 Mpro: an in silico approach towards drug repurposing.

SARS-CoV-2 Mpro is one of the most vital enzymes of the new coronavirus-2 (SARS-CoV-2) and is a crucial target for drug discovery. Unfortunately, there is not any potential drugs available to combat the action of SARS-CoV-2 Mpro. Based on the reports HIV-protease inhibitors can be applied against the SARS by targeting the SARS-CoV-1 Mpro, we have chosen few clinically trialed experimental and allophenylnorstatine (APNS) containing HIV-protease inhibitors (JE-2147, JE-533, KNI-227, KNI-272 & KNI-1931), to examine their binding affinities with SARS-CoV-2 Mpro and to assess their potential to check for a possible drug candidate against the protease. Here, we have chosen a methodology to understand the binding mechanism of these five inhibitors to SARS-CoV-2 Mpro by merging molecular docking, molecular dynamics (MD) simulation and MM-PBSA based free energy calculations. Our estimations disclose that JE-2147 is highly effective (ΔGBind = -28.31 kcal/mol) due to an increased favorable van der Waals (ΔEvdw) interactions and decreased solvation (ΔGsolv) energies between the inhibitor and viral protease. JE-2147 shows a higher level of interactions as compared to JE-533 (-6.85 kcal/mol), KNI-227 (-18.36 kcal/mol), KNI-272 (-15.69 kcal/mol) and KNI-1931 (-21.59 kcal/mol) against SARS-CoV-2 Mpro. Binding contributions of important residues (His41, Met49, Cys145, His164, Met165, Glu166, Pro168, Gln189, etc.) from the active site or near the active site regions with ≥1.0 kcal/mol suggest a potent binding of the inhibitors. It is anticipated that the current study of binding interactions of these APNS containing inhibitors can pitch some valuable insights to design the significantly effective anti-SARS-CoV-2 Mpro drugs.Communicated by Ramaswamy H. Sarma.

In Silico Optimization of SARS-CoV-2 Spike Specific Nanobodies.

BACKGROUND: The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide, caused a global pandemic, and killed millions of people. The spike protein embedded in the viral membrane is essential for recognizing human receptors and invading host cells. Many nanobodies have been designed to block the interaction between spike and other proteins. However, the constantly emerging viral variants limit the effectiveness of these therapeutic nanobodies. Therefore, it is necessary to find a prospective antibody designing and optimization approach to deal with existing or future viral variants. METHODS: We attempted to optimize nanobody sequences based on the understanding of molecular details by using computational approaches. First, we employed a coarse-grained (CG) model to learn the energetic mechanism of the spike protein activation. Next, we analyzed the binding modes of several representative nanobodies with the spike protein and identified the key residues on their interfaces. Then, we performed saturated mutagenesis of these key residue sites and employed the CG model to calculate the binding energies. RESULTS: Based on analysis of the folding energy of the angiotensin-converting enzyme 2 (ACE2) -spike complex, we constructed a detailed free energy profile of the activation process of the spike protein which provided a clear mechanistic explanation. In addition, by analyzing the results of binding free energy changes following mutations, we determined how the mutations can improve the complementarity with the nanobodies on spike protein. Then we chose 7KSG nanobody as a template for further optimization and designed four potent nanobodies. Finally, based on the results of the single-site saturated mutagenesis in complementarity determining regions (CDRs), combinations of mutations were performed. We designed four novel, potent nanobodies, all exhibiting higher binding affinity to the spike protein than the original ones. CONCLUSIONS: These results provide a molecular basis for the interactions between spike protein and antibodies and promote the development of new specific neutralizing nanobodies.

Interfacial water molecules contribute to antibody binding to the receptor-binding domain of SARS-CoV-2 spike protein.

Antibodies that recognize the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), especially the neutralizing antibodies, carry great hope in the treatment and final elimination of COVID-19. Driven by a synchronized global effort, thousands of antibodies against the spike protein have been identified during the past two years, with the structural information available at atomistic detail for hundreds of these antibodies. We developed an improved molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method including explicitly treated interfacial water to calculate the binding free energy between representative antibodies and the receptor binding domain (RBD) domain of SARS-COV-2 spike proteins. We discovered that explicit treatment of water molecules located at the interface between RBD and antibody effectively improves the results for the WT and variants of concern (VOC) systems. Interfacial water molecules, together with surface and internal water molecules, behave drastically from bulk water and exert peculiar impacts on protein dynamics and energy, and thus warrant explicit treatment to complement implicit solvent models. Our results illustrate the importance of including interfacial water molecules to approach efficient and reliable prediction of binding free energy.Communicated by Ramaswamy H. Sarma.

Vitamin D3 and its hydroxyderivatives as promising drugs against COVID-19: a computational study.

The epidemiologic correlation between the poor prognosis of SARS-CoV-2 infection and vitamin D deficiency has been observed worldwide, however, their molecular mechanisms are not fully understood. In this study, we used combined molecular docking, molecular dynamics simulations and binding free energy analyses to investigate the potentials of vitamin D3 and its hydroxyderivatives as TMPRSS2 inhibitor and to inhibit the SARS-CoV-2 receptor binding domain (RBD) binding to angiotensin-converting enzyme 2 (ACE2), as well as to unveil molecular and structural basis of 1,25(OH)2D3 capability to inhibit ACE2 and SARS-CoV-2 RBD interactions. The results show that vitamin D3 and its hydroxyderivatives are favorable to bind active site of TMPRSS2 and the binding site(s) between ACE2 and SARS-CoV2-RBD, which indicate that vitamin D3 and its biologically active hydroxyderivatives can serve as TMPRSS2 inhibitor and can inhibit ACE2 binding of SARS-CoV-2 RBD to prevent SARS-CoV-2 entry. Interaction of 1,25(OH)2D3 with SARS-CoV-2 RBD and ACE2 resulted in the conformation and dynamical motion changes of the binding surfaces between SARS-CoV-2 RBD and ACE2 to interrupt the binding of SARS-CoV-2 RBD with ACE2. The interaction of 1,25(OH)2D3 with TMPRSS2 also caused the conformational and dynamical motion changes of TMPRSS2, which could affect TMPRSS2 to prime SARS-CoV-2 spike proteins. Our results propose that vitamin D3 and its biologically active hydroxyderivatives are promising drugs or adjuvants in the treatment of COVID-19. Communicated by Ramaswamy H. Sarma.

Design, synthesis, and molecular dynamics simulation studies of quinoline derivatives as protease inhibitors against SARS-CoV-2.

A new series of quinoline derivatives has been designed and synthesized as probable protease inhibitors (PIs) against severe acute respiratory syndrome coronavirus 2. In silico studies using DS v20.1.0.19295 software have shown that these compounds behaved as PIs while interacting at the allosteric site of target Mpro enzyme (6LU7). The designed compounds have shown promising docking results, which revealed that all compounds formed hydrogen bonds with His41, His164, Glu166, Tyr54, Asp187, and showed π-interaction with His41, the highly conserved amino acids in the target protein. Toxicity Prediction by Komputer Assisted Technology results confirmed that the compounds were found to be less toxic than the reference drug. Further, molecular dynamics simulations were performed on compound 5 and remdesivir with protease enzyme. Analysis of conformational stability, residue flexibility, compactness, hydrogen bonding, solvent accessible surface area (SASA), and binding free energy revealed comparable stability of protease:5 complex to the protease: remdesivir complex. The result of hydrogen bonding showed a large number of intermolecular hydrogen bonds formed between protein residues (Glu166 and Gln189) and ligand 5, indicating strong interaction, which validated the docking result. Further, compactness analysis, SASA and interactions like hydrogen-bonding demonstrated inhibitory properties of compound 5 similar to the existing reference drug. Thus, the designed compound 5 might act as a potential inhibitor against the protease enzyme.Communicated by Ramaswamy H. SarmaHighlightsQuinoline derivatives have been designed as protease inhibitors against SARS-CoV-2.The compounds were docked at the allosteric site of SARS-CoV-2-Mpro enzyme (PDB ID: 6LU7) to study the stability of protein-ligand complex.Docking studies indicated the stable ligand-protein complexes for all designed compounds.The Toxicity Prediction by Komputer Assisted Technology protocol in DS v20.1.0.19295 software was used to evaluate the toxicity of the designed quinoline derivatives.Molecular dynamics studies indicated the formation of stable ligand-Mpro complexes.

Computational identification of potential inhibitory compounds in Indian medicinal and aromatic plant species against major pathogenicity determinants of SARS-CoV-2.

SARS-CoV-2 (COVID-19) viral pandemic has been reported across 223 countries and territories. Globalized vaccination programs alongside administration of repurposed drugs will assumingly confer a stronger and longer individual specific immune protection. However, considering possible recurrence of the disease via new variants, a conveniently deliverable phytopharmaceutical drug might be the best option for COVID-19 treatment. In the current study, the efforts have been made to identify potential leads for inhalation therapy as nasal swabs have been reported to transfer viral load prominently. In that direction, 2363 Essential oil (EOs) compounds from Indian medicinal and aromatic plants were screened through docking analysis and potential candidates were shortlisted that can interfere with viral pathogenicity. The main protease (Mpro) of SARS-CoV-2 interacted closely with jatamansin (JM), 6,7-dehydroferruginol (FG) and beta-sitosterol (BS), while Papain-like Protease (PLpro) with friedelane-3-one (F3O) and lantadene D (LD) independently. Reduced Lantadene A (LAR) exhibited preferable interaction with RNA-dependent-RNA-polymerase (RdRp) whereas Lantadene A (LA) with RdRp and spike-glycoprotein (SG-pro) both target proteins. When compared against highest binding affinity conformations of well-known inhibitors of targets, these prioritized compounds conferred superior or comparable SARS-CoV-2 protein inhibition. Additionally, promising results were noted from pharmacokinetics prediction for all shortlisted compounds. Besides, molecular dynamics simulation for 100 ns in two replicates and binding free energy analysis revealed the stability of complexes with optimum compactness. To the best of our knowledge, the current investigation is a unique initial attempt whereby EO compounds have been computationally screened, irrespective of their known medicinal properties to fight COVID-19 infection.Communicated by Ramaswamy H. Sarma.

Finding a prospective dual-target drug for the treatment of coronavirus disease by theoretical study.

Spike protein of coronavirus is a key protein in binding and entrance of virus to the human cell via binding to the receptor-binding domain (RBD) domain of S1 subunit to peptidase domain region of ACE2 receptor. In this study, the possible effect of 24 antiviral drugs on the RBD domain of spike protein was investigated via docking and molecular dynamics simulation for finding a dual-target drug. At first, all drugs were docked to the RBD domain of spike protein, and then all complexes and free RBD domains were separately used for molecular dynamics simulation for 50 ns via amber18 software. The simulation results showed that 10 ligands from 28 ligands were separated from the RBD domain, and among 18 remained ligands, baloxavir marboxil, and danoprevir drugs, besides endonuclease activity and protease inhibitory, can bind to key residues of the RBD domain. Then these drugs have a dual target and should be more effective than current drugs, and experimental studies should be done on baloxavir marboxil and danoprevir as more potential drugs for coronavirus disease Communicated by Ramaswamy H. Sarma.

Virtual screening and structure-based 3D pharmacophore approach to identify small-molecule inhibitors of SARS-CoV-2 Mpro.

The outbreak caused by a coronavirus 2 has required quick and potential treatment strategies. The main protease enzyme Mpro plays an important role in the viral replication which renders it an important target for discovering SARS-CoV-2 inhibitors. In this study, 3D pharmacophore structure-based virtual screening and molecular docking were conducted using MOE and Bristol University Docking Engine (BUDE). Around 400,000 molecules of ZINC15 database were docked against the crystal structure of main protease, followed by 3D pharmacophore filtration. Six top-ranked hits (ZINC58717986, ZINC60399606, ZINC58662884, ZINC45988635, ZINC54706757 and ZINC17320595) were identified based on their strong spatial affinity and forming H-bonds with key residues H41, E166, Q189 and T190 of the binding pocket of Mpro SARS-CoV-2. The 6 hits subjected to molecular dynamics simulations for 100 ns followed by binding free energy calculations using MM-PBSA technique. Interestingly, three hits showed free binding energy (ΔGbinding) lower than tert-butyl N-[1-[(2S)-1-[[(2S)-4-(benzylamino)-3,4-dioxo-1-[(3S)-2-oxopyrrolidin-3-yl]butan-2-yl]amino]-3-cyclopropyl-1-oxopropan-2-yl]-2-oxopyridin-3-yl]carbamate (α-ketoamide 13 b) (ΔGbinding) -76.67 ± 0.5 kJ/mol which suggested their potential against SARS-CoV-2. The best binding free energy candidates, ZINC58717986 (ΔGbinding) -98.41 ± 0.7 kJ/mol. The second-best hit candidate, ZINC54706757 (ΔGbinding) -83.4 ± 0.6 kJ/mol and the third one ZINC17320595 (ΔGbinding) -78.85 ± 0.5 kJ/mol. Per residue decomposition free energy indicates H41, S46, H164, E166, D187, Q189 and Q192 are hot spot residues while residues M49, M165, L167 and P168 contribute to the hydrophobic interactions. The pharmacokinetic study suggests that the selected 6 hits possess drug-like properties. The 3D pharmacophore virtual screening, molecular dynamics and MM-PBSA approaches facilitated identification 3 promising hits with low free binding energy as SARS-CoV-2 inhibitors.Communicated by Ramaswamy H. Sarma.

Elucidation of the inhibitory activity of plant-derived SARS-CoV inhibitors and their potential as SARS-CoV-2 inhibitors.

Several drugs are now being tested as possible therapies due to the necessity of treating SARS-CoV-2 infection. Although approved vaccines bring much hope, a vaccination program covering the entire global population will take a very long time, making the development of effective antiviral drugs an effective solution for the immediate treatment of this dangerous infection. Previous studies found that three natural compounds, namely, tannic acid, 3-isotheaflavin-3-gallate and theaflavin-3,3-digallate, are effective proteinase (3CLpro) inhibitors of SARS-CoV (IC50 <10 µM). Based on this information and due to the high sequence identity between SARS-CoV and SARS-CoV-2 3CLpro, these three compounds could be candidate inhibitors of SARS-CoV-2 3CLpro. This paper explores the structural and energetic features that guided the molecular recognition of these three compounds for dimeric SARS-CoV-2 and SARS-CoV 3CLpro, the functional state of this enzyme, using docking and MD simulations with the molecular mechanics-generalized-born surface area (MMGBSA) approach. Energetic analysis demonstrated that the three compounds reached good affinities in both systems in the following order: tannic acid > 3-isotheaflavin-3-gallate > theaflavin-3,3-digallate. This tendency is in line with that experimentally reported between these ligands and SARS-CoV 3CLpro. Therefore, tannic acid may have clinical usefulness against COVID-19 by acting as a potent inhibitor of SARS-CoV-2 3CLpro.Communicated by Ramaswamy H. Sarma.

Multidimensional virtual screening approaches combined with drug repurposing to identify potential covalent inhibitors of SARS-CoV-2 3CL protease.

The outbreak of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has caused an unprecedented global pandemic, and new cases are still on the rise due to the absence of effective medicines. However, developing new drugs within a short time is extremely difficult. Repurposing the existing drugs provides a fast and effective strategy to identify promising inhibitors. Here we focus on the screening of drugs database for discovering potential covalent inhibitors that target 3-chymotrypsin-like protease (3CLpro), an essential enzyme mediating viral replication and transcription. Firstly, we constructed a receptor-ligand pharmacophore model and verified it through decoy set. The importance of pharmacophore features was evaluated by combining molecular dynamics simulation with interaction analyses. Then, covalent docking was used to perform further screening. According to docking score and Prime/Molecular Mechanics Generalized Born Surface Area (MM-GBSA) score, total ten compounds obtained good scores and successfully established covalent bonds with the catalytic Cys145 residue. They also formed favorable interactions with key residues in active sites and closely integrated with 3CLpro with binding modes similar to known 3CLpro inhibitor. Finally, the top four hits DB08732, DB04653, DB01871 and DB07299 were further subjected to 100 ns molecular dynamics (MD) simulation and MM-GBSA binding free energy calculations. The results suggest that the four candidates show good binding affinities for 3CLpro, which warrants further evaluation for their in-vitro/in-vivo activities. Overall, our research methods provide a valuable reference for discovering promising inhibitors against SARS-CoV-2 and help to fight against the epidemic.Communicated by Ramaswamy H. Sarma.

Unveiling the Inhibitory Potentials of Peptidomimetic Azanitriles and Pyridyl Esters towards SARS-CoV-2 Main Protease: A Molecular Modelling Investigation.

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is responsible for COVID-19, which was declared a global pandemic in March 2020 by the World Health Organization (WHO). Since SARS-CoV-2 main protease plays an essential role in the virus's life cycle, the design of small drug molecules with lower molecular weight has been a promising development targeting its inhibition. Herein, we evaluated the novel peptidomimetic azatripeptide and azatetrapeptide nitriles against SARS-CoV-2 main protease. We employed molecular dynamics (MD) simulations to elucidate the selected compounds' binding free energy profiles against SARS-CoV-2 and further unveil the residues responsible for the drug-binding properties. Compound 8 exhibited the highest binding free energy of -49.37 ± 0.15 kcal/mol, followed by compound 7 (-39.83 ± 0.19 kcal/mol), while compound 17 showed the lowest binding free energy (-23.54 ± 0.19 kcal/mol). In addition, the absorption, distribution, metabolism, and excretion (ADME) assessment was performed and revealed that only compound 17 met the drug-likeness parameters and exhibited high pharmacokinetics to inhibit CYP1A2, CYP2C19, and CYP2C9 with better absorption potential and blood-brain barrier permeability (BBB) index. The additional intermolecular evaluations suggested compound 8 as a promising drug candidate for inhibiting SARS-CoV-2 Mpro. The substitution of isopropane in compound 7 with an aromatic benzene ring in compound 8 significantly enhanced the drug's ability to bind better at the active site of the SARS-CoV-2 Mpro.

Bacillus species; a potential source of anti-SARS-CoV-2 main protease inhibitors.

The COVID-19 being a preconized global pandemic by the World Health Organization needs persuasive immediate research for possible medications. The present study was carried out with a specific aim to computationally evaluate and identify compounds derived from Bacillus species as the plausible inhibitors against 3-chymotrypsin-like main protease (3CLpro) or main protease (MPro), which is a key enzyme in the life-cycle of coronavirus. The compounds were isolated from the crude extracts of Bacillus species. Among the isolated compounds, novel inhibitory leads were identified using in silico techniques. Molecular docking revealed that stigmasterol (-8.3 kcal/mol), chondrillasterol (-7.9 kcal/mol) and hexadecnoic acid (-6.9 kcal/mol)) among others bind in the substrate-binding pocket and also interacted with the catalytic dyad of the 3-CLpro. Further evaluation using 50 ns molecular dynamic simulation and MMPB-GBSA indicated that among the top three docking hits, hexadecanoic acid was found to be the most promising anti-COVID-19 lead against the main protease. Hexadecanoic acid might serve as a potent anti-SARS-CoV-2 compound to combat COVID-19, however, in vitro and in vivo validation and optimization is needed.Communicated by Ramaswamy H. Sarma.

Free energy perturbation–based large-scale virtual screening for effective drug discovery against COVID-19

As a theoretically rigorous and accurate method, FEP-ABFE (Free Energy Perturbation-Absolute Binding Free Energy) calculations showed great potential in drug discovery, but its practical application was difficult due to high computational cost. To rapidly discover antiviral drugs targeting SARS-CoV-2 Mpro and TMPRSS2, we performed FEP-ABFE–based virtual screening for ∼12,000 protein-ligand binding systems on a new generation of Tianhe supercomputer. A task management tool was specifically developed for automating the whole process involving more than 500,000 MD tasks. In further experimental validation, 50 out of 98 tested compounds showed significant inhibitory activity towards Mpro, and one representative inhibitor, dipyridamole, showed remarkable outcomes in subsequent clinical trials. This work not only demonstrates the potential of FEP-ABFE in drug discovery but also provides an excellent starting point for further development of anti-SARS-CoV-2 drugs. Besides, ∼500 TB of data generated in this work will also accelerate the further development of FEP-related methods. © The Author(s) 2022.

An insight into SARS-CoV-2 membrane protein interaction with spike, envelope, and nucleocapsid proteins.

Intraviral protein-protein interactions are crucial for replication, pathogenicity, and viral assembly. Among these, virus assembly is a critical step as it regulates the arrangements of viral structural proteins and helps in the encapsulation of genomic material. SARS-CoV-2 structural proteins play an essential role in the self-rearrangement, RNA encapsulation, and mature virus particle formation. In SARS-CoV, the membrane protein interacts with the envelope and spike protein in Endoplasmic Reticulum Golgi Intermediate Complex (ERGIC) to form an assembly in the lipid bilayer, followed by membrane-ribonucleoprotein (nucleocapsid) interaction. In this study, we tried to understand the interaction of membrane protein's interaction with envelope, spike, and nucleocapsid proteins using protein-protein docking. Further, simulation studies were performed up to 100 ns to examine the stability of protein-protein complexes of Membrane-Envelope, Membrane-Spike, and Membrane-Nucleocapsid proteins. Prime MM-GBSA showed high binding energy calculations for the simulated structures than the docked complex. The interactions identified in our study will be of great importance, as it provides valuable insight into the protein-protein complex, which could be the potential drug targets for future studies.Communicated by Ramaswamy H. Sarma.

Use of simplified models for theoretical prediction of the interactions between available antibodies and the receptor-binding domain of SARS-CoV-2 spike protein.

The negative impact of infectious diseases like COVID-19 on public health and the global economy is evident. This pandemic represents a significant challenge for the scientific community to develop new practical analytical methods for accurately diagnosing emerging cases. Due to their selectivity and sensitivity, new methodologies based on antigen/antibody interactions to detect COVID-19 biomarkers are necessary. In this context, the theoretical, computational modeling reduces experimental efforts and saves resources for rational biosensor design. This study proposes using molecular dynamics to predict the interactions between the Receptor Binding Domain (RBD) of the SARS-CoV-2 spike protein simplified model and a set of highly characterized antibodies. The binding free energy of the antigen/antibody complexes was calculated for the simplified models and compared against the complete SARS-CoV-2 ectodomain to validate the methodology. The structural data derived from our molecular dynamics and end-point free energy calculations showed a positive correlation between both approximations, with a 0.82 Pearson correlation coefficient; t = 3.661, df = 3, p-value = 0.03522, with a 95% confident interval. Furthermore, we identified the interfacial residues that could generate covalent bonds with a specific chemical surface without perturbing the binding dynamics to develop highly sensitive and specific diagnostic devices. Communicated by Ramaswamy H. Sarma.

Implications of SARS-CoV-2 spike protein interactions with Zn-bound form of ACE2: a computational structural study.

The COVID-19 pandemic has generated a major interest in designing inhibitors to prevent SARS-CoV-2 binding on host cells to protect against infection. One promising approach to such research utilizes molecular dynamics simulation to identify potential inhibitors that can prevent the interaction between spike (S) protein on the virus and angiotensin converting enzyme 2 (ACE2) receptor on the host cells. In these studies, many groups have chosen to exclude the ACE2-bound zinc (Zn) ion, which is critical for its enzymatic activity. While the relatively distant location of Zn ion from the S protein binding site (S1 domain), combined with the difficulties in modeling this ion has motivated the decision of exclusion, Zn can potentially contribute to the structural stability of the entire protein, and thus, may have implications on S protein-ACE2 interaction. In this study, the authors model both the ACE2-S1 and ACE2-inhibitor (mAb) system to investigate if there are variations in structure and the readouts due to the presence of Zn ion. Although distant from the S1 or inhibitor binding region, inclusion/exclusion of Zn has statistically significant effects on the structural stability and binding free energy in these systems. In particular, the binding free energy of the ACE2-S1 and ACE2-inhibitor structures is - 3.26 and - 14.8 kcal/mol stronger, respectively, in the Zn-bound structure than in the Zn-free structures. This finding suggests that including Zn may be important in screening potentially inhibitors and may be particularly important in modeling monoclonal antibodies, which may be more sensitive to changes in antigen structure.

Computational approaches to define poncirin from Magnolia champaka leaves as a novel multi-target inhibitor of SARS-CoV-2.

Phytochemical-based drug discovery against the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has been the focus of the current scenario. In this context, we aimed to perform the phytochemical profiling of Magnolia champaka, an evergreen tree from the Magnoliaceae family, in order to perform a virtual screening of its phytoconstituents against different biological targets of SARS-CoV-2. The phytochemicals identified from the ethanol extract of M. champaka leaves using liquid chromatography-mass spectroscopy (LC-MS) technique were screened against SARS-CoV-2 spike glycoprotein (PDB ID: 6M0J), main protease/Mpro (PDB ID: 6LU7), and papain-like protease/PLpro (PDB ID: 7CMD) through computational tools. The experimentation design included molecular docking simulation, molecular dynamics simulation, and binding free energy calculations. Through molecular docking simulation, we identified poncirin as a common potential inhibitor of all the above-mentioned target proteins. In addition, molecular dynamics simulations, binding free energy calculations, and PCA analysis also supported the outcomes of the virtual screening. By the virtue of all the in silico results obtained, poncirin could be taken for in vitro and in vivo studies in near future.Communicated by Ramaswamy H. Sarma.